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Abstract 

The  effect  of  scintillation,  arising  from  propagation  through  atmospheric  turbulence,  on  the  sift 
and  error  probabilities  of  a  quantum  key  distribution  (QKD)  system  that  uses  the  weak  laser  pulse 
version  of  the  Bennett-Brassard  1984  (BB84)  protocol  is  evaluated.  Two  earth-space  scenarios 
are  examined:  satellite-to-ground  and  ground- to-satellite  transmission.  Both  lie  in  the  far-field 
power  transfer  regime.  This  work  complements  previous  analysis  of  turbulence  effects  in  near- 
field  terrestrial  BB84  QKD  [Phys.  Rev.  A  67,  022309  (2003)].  More  importantly,  it  shows  that 
scintillation  has  virtually  no  impact  on  the  sift  and  error  probabilities  in  earth-space  BB84  QKD, 
something  that  has  been  implicitly  assumed  in  prior  analyses  for  that  application.  This  result 
contrasts  rather  sharply  with  what  is  known  for  high-speed  laser  communications  over  such  paths, 
in  which  deep,  long-lived  scintillation  fades  present  a  major  challenge  to  high-reliability  operation. 
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I.  INTRODUCTION 


In  [1]  we  used  the  normal-inode  decomposition  for  line-of-sight  optical  propagation 
through  atmospheric  turbulence  [2]  to  obtain  upper  and  lower  bounds  on  the  sift  and  er¬ 
ror  probabilities  of  a  free-space  optical  implementation  [3]  of  the  Bennett-Brassard  1984 
(BB84)  protocol  for  quantum  key  distribution  (QKD).  There  we  focused  on  modest-length 
terrestrial  paths,  for  which  the  transmit  and  receive  pupils  could  easily  be  large  enough 
that  operation  is  in  the  near-field  power  transfer  regime.  Earth-space  BB84  QKD  is  being 
considered  for  transcontinental  to  worldwide  applications,  and  these  configurations  will  be 
in  the  far-field  power  transfer  regime.  To  date,  almost  all  assessments  of  the  impact  of  at¬ 
mospheric  propagation  on  such  far-field  systems  have  not  included  scintillation  [4],  i.e. ,  only 
the  effects  of  atmospheric  extinction,  turbulence-induced  beam  spread  and  angular  spread, 
and  background-light  collection  have  been  quantified,  see,  e.g.,  [5].  In  this  paper  we  will 
rectify  that  deficiency.  In  particular,  we  will  show  that  scintillation  has  virtually  no  effect 
on  the  sift  and  error  probabilities  of  a  BB84  QKD  system  that  uses  satellite-to-ground  or 
grouncl-to-satellitc  transmission.  This  result  contrasts  rather  sharply  with  what  is  known 
for  high-speed  laser  communications  over  such  paths,  in  which  deep  and  long-lived  scintil¬ 
lation  fades  present  a  major  challenge  arguably  the  major  challenge  to  high-reliability 
operation. 

The  rest  of  this  paper  is  organized  as  follows.  In  Sec.  II  we  briefly  describe  the  BB84 
QKD  protocol  that  we  will  treat,  and  give  expressions  for  its  sift  and  error  probabilities 
conditioned  on  knowledge  of  the  fractional  transmitter-to-receiver  power  transfer.  These 
conditional  probabilities  are  independent  of  whether  operation  is  in  the  near-field  or  far-field 
power  transfer  regimes,  so  they  are  available  from  [1]  for  the  weak  laser  pulse  system  that  we 
shall  consider.  In  Sec.  Ill  we  address  the  satellite- to-ground  link.  Using  a  lognormal-fading 
model  we  show  that  aperture-averaged  scintillation  has  negligible  effect  on  the  uncondi¬ 
tional  sift  and  error  probabilities  [6],  Section  IV  presents  a  similar  analysis  for  two  versions 
of  the  ground-to-satellite  link.  The  first  employs  a  collimated-beam  (nonaclaptive)  transmit¬ 
ter.  The  second  is  an  ideal  adaptive-optics  transmitter  that  perfectly  tracks  the  maximum 
power-transfer  input  eigenfunction  for  the  ground-to-space  path.  In  both  cases  the  conclu¬ 
sion  reached  is  that  same  as  found  in  the  satellitc-to-ground  analysis,  viz.,  scintillation  has 
negligible  effect  on  the  unconditional  sift  and  error  probabilities. 
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II.  SIFT  AND  ERROR  PROBABILITIES  FOR  BB84  FREE-SPACE  QKD 

The  QKD  system  we  consider  is  the  one  described  in  [1].  It  uses  a  line-of-sight  optical 
link  to  connect  a  transmitter  (Alice,  shown  in  Fig.  1)  with  a  receiver  (Bob,  shown  in  Fig.  2). 
On  each  bit  interval,  Alice  chooses  randomly  between  two  linear  polarization  bases,  0°/90° 
and  =f45°,  which  we  will  denote  +  and  X,  respectively.  Having  chosen  a  basis,  she  sends  a 
random  bit  value,  0  or  1,  using  the  coding, 


0°,  if  +  was  chosen, 
—45°,  if  X  was  chosen, 


(la) 


(lb) 


Bob’s  receiver  uses  a  passive  50-50  beam  splitter  to  create  inputs  for  a  pair  of  polarization 
analysis  systems — one  for  the  +  basis  and  one  for  the  X  basis — that  employ  identical  single- 
photon  avalanche  photodiodes  (APDs),  each  of  quantum  efficiency  ?/  [7].  For  a  single  photon 
arriving  at  Bob’s  receiver,  this  passive  arrangement  amounts  to  a  random  choice  between 
the  +  and  the  X  measurement  bases. 

Let  { N0o ,  N$0o ,  AT450 ,  yV+45o }  denote  the  photon  counts  from  the  four  APDs  during  a 
single  bit  interval.  Bob  has  a  detection  event  when  N0o  +  Ar90°  +  AT45o  +  AT+45o  =  1,  i.e., 
when  exactly  one  of  his  detectors  registers  a  count.  In  the  BB84  protocol,  Bob  discloses 
to  Alice  the  sequence  of  bit  intervals  and  associated  measurement  bases  for  which  he  has 
detections.  Alice  then  informs  Bob  which  detections  occurred  in  bases  coincident  with  the 
ones  that  she  used.  These  are  the  sift  events,  i.e.,  bit  intervals  in  which  Bob  has  a  detection 
and  his  count  has  occurred  in  the  same  basis  that  Alice  used.  For  example,  if  Alice  sent 
her  bit  value  as  a  90°-polarized  laser  pulse,  then  a  sift  event  means  that  Bob  had  detected 
exactly  one  count  from  his  four  detectors,  with  N0»  =  1  or  A7)0°  =  1.  An  error  event  is  a  sift 
event  in  which  Bob  decodes  the  incorrect  bit  value.  For  example,  if  Alice  sent  her  bit  value  as 
a  90°-polarized  laser  pulse,  then  an  error  event  means  that  Bob  had  a  sift  in  which  No °  =  1 
occurred.  Once  sift  events  have  been  identified,  the  remainder  of  the  BB84  protocol — -which 
shall  not  concern  us  in  this  paper  -is  standard.  Alice  and  Bob  follow  a  prescribed  set  of 
operations  to  identify  errors  in  their  sifted  bits,  correct  these  errors,  and  apply  sufficient 
privacy  amplification  to  deny  useful  key  information  to  any  potential  eavesdropper  (Eve). 
At  the  end  of  the  full  QKD  procedure,  Alice  and  Bob  have  a  shared  one-time  pad  with 
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which  they  can  communicate  in  complete  security.  For  a  given  level  of  privacy  amplification 
(security),  the  principal  figure-of- merit  for  the  BB84  QKD  system  is  its  key  rate,  i.e. ,  the 
number  of  one-time  pad  bits/sec  that  Alice  and  Bob  produce.  Key  rate  decreases  with 
decreasing  sift  probability  and  increasing  error  probability.  Our  objective  is  to  determine 
the  degree  to  which  turbulence  affects  these  probabilities  [8]. 


FIG.  1:  Block  diagram  of  a  single-laser  QKD  transmitter  (Alice).  The  laser  output  is  a  stream 
of  linearly-polarized  pulses.  The  polarization  controllers  are  driven  by  a  pair  of  random  binary 
sequences.  The  first  sequence  determines  the  sequence  of  polarization  bases  that  will  be  sent:  + 
=  0°/90°  or  x  =  q:45°.  The  second  sequence  determines  the  bit  value  to  be  sent,  according  to  the 
coding  rule  given  in  Eq.  (1).  The  attenuator  reduces  the  transmitter’s  output  to  ns  photons,  on 
average,  per  bit  interval. 

As  in  [1],  we  shall  assume  that  Alice  transmits  an  appropriately  polarized  laser  signal 
pulse  with  an  average  photon  number  of  ns  to  represent  her  bit  value.  Bob’s  receiver  will 
collect  a  random  fraction,  7,  of  the  transmitted  photons  owing  to  the  combined  effects 
of  diffraction,  atmospheric  turbulence,  and  (absorption-plus-scattering  induced)  extinction. 
Indeed,  because  Bob’s  receiver  will  employ  a  narrow  field  of  view — to  minimize  background 
light  shot  noise — it  will  collect  only  the  turbulence-modified  extinguished  direct  beam  from 
Alice’s  transmitter,  i.e.,  no  scattered  light  will  be  collected.  Moreover,  for  bit  durations 
that  are  appreciably  shorter  than  1  ms  and  appreciably  longer  than  1  ps,  we  can  neglect 
time-dependent  fading  and  multipath  spread,  and,  because  atmospheric  turbulence  is  non¬ 
depolarizing,  we  then  have  that  attenuation  by  the  capture  fraction  7  is  the  only  propagation 
effect  incurred  by  Alice’s  transmitted  pulse  cn  route  to  Bob’s  receiver..  I11  addition,  Bob’s 
receiver  will  collect  ns  background  photons  per  polarization,  on  average,  and  each  of  his 
detectors  will  be  subject  to  a  dark-current-equivalent  average  photon  number  of  ns- 

Again  following  [1],  we  shall  assume  our  detectors  have  photon- number  resolution  ca¬ 
pability  and  that  the  their  dead  time  and  afterpulsing  can  be  neglected,  so  that  they  are 
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FIG.  2:  Block  diagram  of  a  QKD  receiver  (Boh).  IF:  interference  filter,  provides  spectral  discrimi¬ 
nation  against  background  light.  50/50:  ordinary  beam  splitter,  provides  a  passive,  random  choice 
of  polarization-analysis  basis  (+  or  x)  for  a  single  photon.  HWP:  half-wave  plate,  converts  x  basis 
into  +  basis.  PBS:  polarizing  beam  splitter.  FS:  field  stop,  provides  spatial- mode  discrimination 
against  background  light.  APD:  single-photon  (Geiger  mode)  avalanche  photodiode. 

governed  by  conditionally-Poisson  counting  statistics.  We  then  have  that  the  conditional 
sift  and  error  probabilities,  given  the  fractional  power  transfer  7  [9],  are 

Pr( sift  |  7)  =  ?/(??.57/2  +  2 nN)e~’t{ns'l+‘lnN\  (2) 

and 

Pr(  error  |  7 )  =  VnNe-’'{nsl+4,'N) .  (3) 

Taking  0  <  C  <  1  to  be  the  atmospheric  extinction  encountered  along  the  propagation  path, 
wc  write  7  =  //£,  where  /t  differs  from  unity  solely  because  of  atmospheric  turbulence. 

Our  interest  is  in  the  unconditional  sift  and  error  probabilities,  i.e., 

Pr(sift)  =  [  d/.ip(p)i](nspC/2  +  2nN)e-^nsl‘c+4nN\  (4) 

J  0 

and 

Pr(error)  =  f  dpp(p)rjn^e~,^nst‘c+iUN\  (5) 

Jo 

where  p(p)  is  the  probability  density  function  (pdf)  for  the  random  variable  //. 

Consider  far- field  propagation  from  the  ground  (z  =  0)  to  a  satellite  (z  =  L),  or  vice 
versa,  in  the  absence  of  turbulence.  In  either  case  a  collimated-beam  transmitter  will  realize 
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2 


(6) 


a  deterministic  fractional  power  transfer  given  by  [10] 


7nt  —  /t  nt£ 


/  nPGPs\ 
V  4A  L  ) 


C  <  1, 


where  Dc  and  Ds  are  the  diameters  of  the  circular  exit/entrance  ground  and  satellite  pupils, 
A  is  the  laser  wavelength,  L  is  the  path  length, 


C  =  exp 


/  dz  a(z)  ] 


(7) 


\  Jo  / 

gives  the  extinction  along  the  path  in  terms  of  the  local  extinction  coefficient,  0(2),  and  the 
subscript  NT  denotes  “no-turbulence”.  In  this  case  we  have  p(/i)  =  5(p  —  pm'),  where  £(•) 
is  the  impulse  function,  so  that 


Pr(sift)NT  =  //(n5/tNT£/2  +  2  nN)e~^nstwC+4nN\ 


(8) 


and 

Pr(error)NT  =  7 1nNe~^nsllti'rC+inN\  (9) 

are  the  no-turbulence  sift  and  error  probabilities. 

To  proceed  further  we  will  consider  satellite-to-ground  and  ground-to-satellite  scenarios, 
so  that  we  can  employ  specific  pdfs  for  //.  and  evaluate  Pr(sift.)  and  Pr (error). 


III.  SATELLITE-TO-GROUND  SCENARIO 


Consider  the  satellite-to-ground  scenario  in  which  the  satellite’s  diameter- As  exit  pupil 
lies  well  within  a  single  turbulence  coherence  area  but  the  ground  terminal’s  diameter- 
Dq  entrance  pupil  comprises  many  turbulence  coherence  areas.  In  this  case  the  extended 
Huygens- Fresnel  principle  leads  to  the  following  expression  for  the  no-extinction  fractional 
power  transfer,  assuming  that  the  ground  detector’s  field  of  view  is  large  enough  to  capture 
all  the  laser  light  reaching  the  ground  terminal’s  entrance  pupil  [11]: 


li  = 


nDs 

4(AL)2 


j\p\<Da/2 


dpe2x{0'p). 


(10) 


Here,  y(0,  p)  is  the  log-amplitude  fluctuation  imposed  on  the  field  received  at  transverse 
coordinate  p  in  the  2  =  0  plane  from  a  point  source  radiating  from  transverse  coordinate  0 
in  the  2  =  L  plane.  Energy  conservation  implies  that 


where  (•)  denotes  ensemble  average,  whence 


2 


</*>  = 


(  kDgDs\ 

V  4A L  ) 


—  /ii\'T- 


(12) 


Equation  (12)  represents  well-known  downlink  behavior,  i.e. ,  there  is  no  beam  spread  due 
to  turbulence  when  the  transmitter  pupil  lies  within  a  single  coherence  area  so  that  aver¬ 
age  power  transfer  is  unaffected  by  turbulence  if  the  receiver’s  field  of  view  is  sufficient  to 
accommodate  any  turbulence-induced  angular  spread  [12]. 

At  this  point,  we  will  assume  that  y(0,  p)  is  Gaussian  distributed,  statistically  homoge¬ 
neous,  and  isotropic.  Thus  it  is  completely  characterized  by  its  mean  mx,  variance  a2,  and 
covariance  function 


Kxx(p)  =  (AX(0,p)AX(0,0)),  (13) 

with  p  =  \p\  and  Ay(0,  p)  =  y(0,  p)  —  mx.  The  energy  conservation  condition  in  Eq.  (12) 
then  implies  that  mx  =  —  o2v  Furthermore,  because  the  sum  of  real- valued  lognormal 
random  variables  is  well  approximated  by  a  real- valued  lognormal  random  variable  [13],  we 
can  write  [14] 

p  =  pNre2u,  (14) 


where  u  is  a  Gaussian  random  variable  with  mean  —  a'2  and  variance  o20  with 


eAa*  -1=  4 


•dg 


*dg  Jo 


dp  p(e 


Ai^x(p)  _ 


1) 


2  r 


x  — 

7T  L 


cos  \p/Dg)  -  {p/DG)\f 1  -  {p/DcY 


(15) 


Equation  (15)  is  an  aperture-averaging  formula,  viz.,  it  amounts  to 


.An}, 


1  = 


e'1<7x  -  1 

A'v(Dg)’ 


(l(i) 


where  Nx(Dg)  3>  1  bj'  assumption  is  the  number  of  log-amplitude  coherence  areas  in  the 
ground  terminal’s  entrance  pupil. 

Unfortunately,  the  lognormal  distribution  for  p  does  not  permit  us  to  obtain  closed-form 
expressions  for  the  sift  and  error  probabilities.  So,  to  demonstrate  that  scintillation  has 
essentially  no  effect  on  these  probabilities,  we  will  take  a  numerical  approach  employing  a 
worst-case  value  for  a2.  In  particular,  we  shall  choose  a2  —  0.5,  which  is  a  conservative 
upper  bound  on  saturated  scintillation,  and  use  NX(DG )  =  1  in  Eq.  (16),  thus  obtaining  a 
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worst-case  value  er2  —  0.5.  Using  the  reasonable  values  ns  =  0.5,  n,v  =  5  x  10  G,  ij  =  0.5, 


7nt  =  10  3,  we  find 


and 


Pr(sift)NT  =  1.299  x  HT4, 


Pr(error)NT  =  2.499  x  10  6. 

Using  a  lognormal  distribution  for  u  with  mu  =  — cr2  =  —0.5  we  then  find  that 

Pr(sift) 


and 


Pr(sift)NT 

Pr(error) 


-  1 


-  1 


=  1.52  x  1(T3, 


=  1.97  x  IQ"7, 


(17) 

(18) 

(19) 

(20) 


Pr(error)NT 

showing  that  scintillation  has  indeed  had  a  negligible  influence  on  the  sift  and  error  proba¬ 
bilities. 

Before  moving  on  to  the  ground-to-satellite  scenario,  it  is  worth  noting  that  security  of 
BB84  QKD  depends  on  Pr(  error  |  sift ) — the  conditional  error  probability  given  that  a  sift 
event  has  occurred — being  sufficiently  small.  This  quantity  is  usually  called  the  quantum 
bit-error  rate  (QBER).  For  the  example  given  above  we  find  that 

Pr(  error  |  sift)NT  =  1-92  x  10-2, 


in  the  absence  of  turbulence,  and  we  get 

Pr(  error  |  sift ) 


-  1 


=  1.53  x  10"3, 


(21) 


(22) 


Pr(  error  |  sift  )nt 

in  the  presence  of  worst-case  scintillation.  Thus  worst-case  scintillation  on  the  satellite-to- 
ground  path  has  virtually  no  effect  on  the  QBER. 


IV.  GROUND-TO-SATELLITE  SCENARIO 

Now  let  us  turn  to  the  ground-to-satellite  scenario.  We  will  continue  to  make  the  ge¬ 
ometric  assumptions  that  we  employed  for  the  satellite-to-ground  case,  namely  that  the 
ground  terminal’s  diameter- Dq  exit  pupil  comprises  many  turbulence  coherence  areas,  and 
that  the  satellite’s  diameter-  Ds  entrance  pupil  lies  well  within  a  single  turbulence  coherence 
area.  Here,  however,  it  turns  out  that  there  are  two  cases  worth  considering:  nonadap- 
tive  versus  adaptive-optics  transmitters.  VVe  will  start  with  the  nonadaptive  case  using  a 
collimated-beam  transmitter. 
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A.  Collimated-Beam  Transmitter 


In  the  absence  of  turbulence,  a  ground-based  collimated-beam  transmitter  whose  normal¬ 
ized  spatial  mode  pattern  is 


&(p)  = 


for  |p|  <  Do/2, 

0,  otherwise, 


(23) 


achieves  7nt  for  its  far-held  power  transfer  to  the  satellite.  In  the  presence  of  turbulence, 
this  same  collimated-beam  transmitter  yields  7  =  p£  for  its  fractional  power  transfer,  where 


\DC\L) 


L 


dp  ex 


\(0,p)+i<l>(0,p) 


P<Oa/2 


(24) 


In  this  expression,  y(0,  p)  and  0(0,  p)  are  the  log-amplitude  and  phase  fluctuations  imposed 
on  the  Held  received  at  transverse  coordinate  0  in  the  2  =  L  plane  from  a  point  source 
at  transverse  coordinate  p  in  the  2  =  0  plane.  Atmospheric  reciprocity  [15]  implies  that 
this  x'(0,  p)  and  0(0,  p)  can  also  be  regarded  as  the  log-amplitude  and  phase  fluctuations 
imposed  on  the  Held  received  at  transverse  coordinate  p  in  the  2  =  0  plane  from  a  point 
source  radiating  from  transverse  coordinate  0  in  the  2  =  L  plane. 

The  average  behavior  of  the  preceding  p  is  well  known,  see,  e.g.,  [11].  For  our  case,  in 
which  the  ground  terminal’s  exit  pupil  comprises  many  turbulence  coherence  areas,  we  can 
use  the  asymptotic  result, 


</*>  = 


/  7 tDtDs 

V  4A  L 


2 


(25) 


where 


Dr  =  3.18po 


(26) 


gives  the  effective  (turbulence-limited)  transmitter  diameter  in  terms  of  the  turbulence  co¬ 
herence  length 


Po  = 


dzCl{z)(\-z/Lf/*\ 


-3/5 


(27) 


with  k  =  2n/\  being  the  wave  number  at  the  laser  wavelength,  and  Cf,(z)  being  the  tur¬ 
bulence  strength  parameter  along  the  propagation  path  from  the  ground  terminal  to  the 
satellite.  Note  that  this  result,  which  was  originally  derived  in  the  weak-perturbation 
(Rytov-approximation)  regime,  is  valid  in  the  strong-perturbation  regime.  More  impor¬ 
tantly,  Dt  “C  Dq,  because  we  have  assumed  that  the  ground  terminal’s  exit  pupil  comprises 
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a  large  number  of  turbulence  coherence  areas.  It  follow  that.  (/t)  <C  pNT  for  the  collimated- 
bcam  (nonadaptive)  ground-to-satellite  transmitter.  This  is  the  well-known  beam  spread 
result,  which  has  driven  earth-space  BB84  QKD  designs  to  prefer  satellite-to-ground  oper¬ 
ation,  wherein  (p)  =  //ktt  prevails  [5],  as  we  saw  in  Sec.  III. 

To  complete  our  evaluation  of  the  sift  and  error  probabilities  for  the  collimated-beam 
transmitter,  we  make  use  of  the  Central  Limit  Theorem  argument  that  implies  //  will  be 
exponentially  distributed  when  the  ground  terminal’s  exit  pupil  contains  a  large  number 
of  turbulence  coherence  areas.  Unlike  the  lognormal  statistics  that  arose  in  Sec.  Ill,  the 
exponential  distribution  leads  to  simple  closed- form  results,  i.e. , 


Pr(sift) 


ye  4,1,1  N  /  ns(fi)C  4 nN  \ 

2  V  (1  +  Vns{li)£)2  1  +  yns(n)c)  ' 


(28) 


and 


Pr (error)  = 


rtnNe 


(29) 


1  + 1 ins(n)C 

In  seeking  to  assess  the  impact  of  scintillation  on  this  nonadaptive- transmitter  QKD  system, 
it  is  unfair  to  compare  these  probabilities  to  the  results  from  Eqs.  (8)  and  (9),  because  those 
sift  and  error  probabilities  assume  that  there  is  no  turbulence,  whereas  turbulence-induced 
beam  spread — which  is  due  primarily  to  the  phase  fluctuations  0(0,  p) — has  a  dramatic 
impact  on  the  average  ground-to-satellite  power  transfer  even  if  there  are  no  log-amplitude 
fluctuations,  viz.,  x(0,  p)  =  0.  So,  for  our  comparison  case,  we  will  consider  a  ground-to- 
satellite  link  in  which  p(p)  =  d(p  —  {//.))  with  (/i)  given  by  the  turbulence-limited  result  from 
Eq.  (25).  This  no-scintillation  case  then  has  the  following  sift  and  error  probabilities: 


Pr(sift)NS  =  ii(ns(i<)£.ri  +  2n«)e-''<"sWC+J”'',), 


(30) 


and 

Pr  (err  or)  Ns  =  i]riNe~^n^l)c+AnN) .  (31) 

Figure  3  shows  plots  of  Pr(sift)NSi  Pi‘(error)Ns  and  Pr(  error  |  sift  )ns  versus  {p)C.  The 
ns,  njv,  and  //  values  assumed  are  the  same  as  were  used  in  the  satellite-to-ground  example 
in  Sec.  Ill,  viz.,  ns  =  0.5,  n,y  =  5  x  10_(i,  and  //  =  0.5.  These  curves  are  indistinguishable 
from  similar  plots  for  Pr(sift),  Pr(error)  and  Pr( error  |  sift)  because,  over  the  range  10-4  < 
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(//,)£  <10  2,  we  have  that 


Pr(sift) 

Pr(sift)NT 


2.48  x  1()-3, 


(32) 


max 


Pr  (error) 
Pr(error)i\T 


3.11  x  10~G, 


(33) 


when  ns  =  0.5,  n n 
scintillation. 


max 


Pr(  error  |  sift ) 


1 


=  2.49  x  10“3, 


(34) 


Pr(  error  |  sift  )nt 
5  x  10-6,  and  //  =  0.5.  Once  again  we  see  virtually  no  effect  from 


FIG.  3:  (Color  online)  No-scintillation  sift  probability,  error  probability,  and  quantum  bit-error 
rate  plotted  versus  the  logarithm  of  the  average  fractional  power  transfer,  (/r)£.  These  curves 
assume  ns  =  0.5,  njv  =  5  x  10-6,  and  ij  =  0.5. 


B.  Optimum  Adaptive-Optics  Transmitter 


For  the  far-field  ground-to-satellite  scenario  we  are  considering,  it  is  known  that  the 
optimum  adaptive-optics  transmitter — which  perfectly  adapts  both  the  amplitude  and  the 
phase  of  the  transmitter’s  spatial  mode  pattern  to  maximize  the  fractional  power  transfer 
produces  the  following  normalized  spatial  mode  pattern  [10], 


€(P) 


,x(Op)-j<P(U-P) 


.  (L 


dp  e  x 


2\(0  .p) 


\p<Da/2 


7721  for  \p\  <  Dc/2 , 


(35) 


0, 


otherwise. 
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This  transmitter  achieves 


h  = 


7tD| 

4(AZ)s 


<\p\<Dc/2 


dp  e  x 


2.V(0  ,p) 


(36) 


because  it  is  the  reciprocity  dual  of  the  satellite-to-ground  case.  Hence  the  same  considera¬ 
tions  made  in  Sec.  Ill  show  that  scintillation  has  essentially  no  effect  on  the  sift  probability, 
error  probability,  and  quantum  bit-error  rate  of  this  system. 


V.  CONCLUSIONS 

We  have  evaluated  the  sift  probability,  error  probability,  and  quantum  bit-error  rate 
for  satellite-to-ground  and  ground-to-satellitc  BB84  QKD  when  the  ground  terminal’s  en¬ 
trance/exit  pupil  contains  a  large  number  of  turbulence  coherence  areas,  the  satellite’s  en¬ 
trance/exit  pupil  lies  well  within  a  single  turbulence  coherence  area,  and  operation  is  deep 
into  the  far-field  power  transfer  regime.  For  reasonable  choices  of  the  average  transmit¬ 
ted  photon  number,  the  average  number  of  noise  photons  reaching  each  detector,  and  the 
quantum  efficiency,  we  have  found  that  scintillation  has  no  appreciable  effect  on  the  afore¬ 
mentioned  performance  metrics.  Two  final  questions  are  worth  addressing.  First,  why  is  it 
that  scintillation  is  so  impotent  here,  when  it  has  long  been  known  to  have  major  impact  on 
the  error  probabilities  of  uncoded  laser  communication  systems?  Second,  would  our  results 
change  were  we  to  use  the  gamma-gamma  distribution,  instead  of  the  lognormal  distribu¬ 
tion,  for  the  satellite-to-ground  scenario  and  the  optimum  adaptive  optics  ground-satellite 
scenario?  The  reason  for  the  disparity  in  scintillation  effects  is  that  laser  communication 
systems  are  trying  to  operate  at  very  low  error  probabilities,  whereas  the  typical  QBER  of 
a  frec-space  optical  BB84  QKD  system  is  a  few  percent.  As  a  result,  it  is  hard  for  even  the 
most  severe  fading  to  change  the  QKD  performance  from  what  is  achieved  at  the  average 
value  of  the  fractional  power  transfer.  Consequently,  we  do  not  expect  the  results  we  re¬ 
ported  for  the  lognormal  distribution  in  Sec.  Ill  to  be  changed  in  any  significant,  way  if  we 
did  a  similar  calculation  using  the  gamma-gamma  distribution. 
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